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The task of analytically diagonalizing a tridiagonal matrix can be considerably simplified when 
a part of the matrix is uniform. Such quasi-uniform matrices occur in several physical contexts, 
both classical and quantum, where one-dimensional interactions prevail. These include magnetic 
chains, ID arrays of Josephson junctions or of quantum dots, boson and fermion hopping models, 
random walks, and so on. In such systems the bulk interactions are uniform, and differences may 
occur around the boundaries of the arrays. Since in the uniform case the spectrum consists of a 
band, we exploit the bulk uniformity of quasi-uniform tridiagonal matrices in order to express the 
spectral problem in terms of a variation of the distribution of eigenvalues in the band and of the 
■ corresponding eigenvectors. In the limit of large matrices this naturally leads to a deformation of 

the density of states which can be expressed analytically; a few out-of-band eigenvalues can show up 
' and have to be accounted for separately. The general procedure is illustrated with some examples. 

(N" 



I. INTRODUCTION 

The diagonalization of quasi-uniform tridiagonal matrices, namely tridiagonal matrices which are uniform except 
at the boundaries, appears in many branches of mathematics and physics ^n^j. In particular, tridiagonal matrices 
generally occur in the theory of one-dimensional lattices with nearest-neighbour interactions. In this context quasi- 
• uniform tridiagonal matrices have been recently applied for achieving high quality quantum communication between 
distant parts 611], and for describing spin systems in a spin environment 12 14]. 
r~| ■ In this paper we put forward a general method for calculating the eigenvalues and the eigenvectors of symmetric 
" tridiagonal matrices by exploiting the property of bulk uniformity. This allows us to put the eigenvalues in the form 

a of deformations, defined by suitable shifts, of those of the fully uniform case, which are known to form a band. The 
' modified density of the eigenmodes in the band is expressed in terms of functions which can be analytically evaluated 

and depend on the non-uniform matrix elements. A small number of localized eigenstates could emerge from the band 
and have to be accounted for separately. Particular applications of this technique were made in Refs. [1, HI]. 

> 

O ■ II. TRIDIAGONAL MATRICES 

m 

A symmetric Ixl tridiagonal matrix T = {T^^} has 2^—1 independent real elements, namely T^^ = (/i = 1, . . . , 
. and Tp^p+i = T^+i^^ = 6^ = 1, . . . , ^— 1). Its spectral decomposition is T = O^^ A O, where O — {Ofc^} is orthogonal, 
■ its rows being the (. eigenvectors of T with eigenvalues Afc, and A = diag({Afc}). 

T is said to be mirror- symmetric if it is also symmetric with respect to the skew diagonal, namely JT, J] = , 
where J^^ ~ 5f^j+i-v is the mirroring matrix. In the mathematical language such matrices are both persymmetric 
(JTJ — T*) and centra symmetric {JTJ = T). It is known that the eigenvectors of a mirror-symmetric T are either 
. symmetric or antisymmetric [l6j . 

^ S : Ofc,£+i-M = ; (1) 

this formula assumes that &^ > and the eigenvalues {Afc} are listed in decreasing order. 

The eigenvectors can be completely expressed in terms of characteristic polynomials of submatrices of T, evaluated 
at the eigenvalues. In order to prove this, let us introduce the following notation for tridiagonal submatrices. 



bfi a^+i 6p+i 



T 



ay J 



(2) 
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and for the corresponding characteristic polynomials, 

X,.,.(A)=dct[A-r^:.] , (3) 

where /i < i^; then Ti-j = T and Xi-.eW = x(A); the eigenvalues are the £ solutions of the secular equation xi^k) — 0. 
By expanding from the bottom (upper) corners, these polynomials are found to satisfy the recurrence relations 



X.:.(A) - (A - a.) x,:.-,(A) - bl_, x,:.-.(A), (4a) 



X,:„ (A) = (A - a^) x,+.. (A) - bl (A) . (4b) 



The following important and useful formula (see, e.g., Ref. [131) expresses the product of two components of the 
same eigenvector, 



X„ + l:f(Afc) , (5) 



which holds for /i < if one defines Xi,o(Afe) = X^+i f (Afc) = 1- One can assume 6^ 7^ for all /j. = 1, 
as otherwise the diagonalization of T would split into the diagonalization of independent submatrices, so that the 
eigenvalues of T are nondegenerate. Hence, the derivatives of the characteristic polynomial at the eigenvalues do not 
vanish, x'(Afe) ^ 0, and Eq. ([5]) can be solved for the eigenvector components, for example 



.2 _ X2:«(Afe) ^2 _ Xl:<!-l(Afc) 



- X'(A.) ' - X'(A,) ' 

and from one of these (one can arbitrarily choose the positive root) the remaining elements of the fc-th eigenvector 
follow by means of Eq. ([S]); for instance, taking = 1, 

r\ r\ I, I, X.+i:t{^k) „ T-r , X^,+i:«(Afe) , , 

Oky ^ Oki bl- ■ -b^-i 7T-^ = Ofci \\b^,-l— — j^-- ■ (7) 

X2:AAfe) f;^^ X^:AAfc) 

This shows that the orthogonal matrix O can be fully expressed in terms of characteristic polynomials. 
Note also that the recurrence equations ^ give 

Ofc..+i = ^^7^0fe.- V^Ofe.,_i (^ = 1,...,^-1) , (8) 

Ou by 

with the assumption Oko — 0; these equations can be used for a sequential computation of the eigenvectors' components 
once the eigenvalues are known. An important consequence of this construction is that, once the first components of 
the eigenvectors {Oki) are determined from Eq. the eigenvectors come out already normalized, i.e., the matrix O is 
orthogonal, making the explicit normalization unnecessary and tremendously simplifying the analytical calculations. 

When the matrix size £ is large the characteristic polynomials Xf,.„(Afc) have a high degree and the analytical 
evaluation of the eigenvalue decomposition is very demanding. In the forthcoming section we provide general simplified 
formulas for the eigenvalues and for the eigenvector elements, Eqs. ^ and ([T]), in the case of a quasi- uniform matrix 
T. 



III. QUASI-UNIFORM TRIDIAGONAL MATRICES 
A. Uniform tridiagonal matrices 

A uniform tridiagonal matrix has equal elements within each diagonal, namely = a and b^ = b, and without loss 
of generality one can set 6 = 1 and a — 0. In this case the recurrence relations (|H) for the characteristic polynomials 
are found to be equal to those defining the Chebyshev polynomials of the second kind [l8l |. 



l^niO , (9) 
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the correspondence being Xi-tW = Ui{\/2). Setting ^ = 
compactly written as 

Un{cosk) - 

so that the secular equation x(A) — Ui{\/2) — defines 

^^^^ = ^' (11) 

With no ambiguity we will use henceforth the index k as running over such a set of £ discrete values, so we may 
keep the notations introduced for the spectral decomposition and, e.g., write the eigenvectors of the uniform case as 
Ofc^ = ^2/{l+l) sMpk). 



cos k the Chebyshev polynomials of the second kind can be 
sin[(n+l)fc] 



sin k ' 

the i eigenvalues A = 2 cos k corresponding to 



(10) 



B. Quasi-uniform tridiagonal matrices 

A tridiagonal matrix T is said to be quasi uniform if it is mainly constituted by a large uniform tridiagonal block 
Tu:y of size nxn (with n = v—u+l), i.e., its elements are a„ = au+i = ■ ■ ■ = = a and 6^ = hu+i = • • • = hy-i = h. 
By 'large uniform block' it is meant that the number of different elements, sitting at one or both corners, is much 
smaller than the size of the whole matrix T, namely that ^—n <C (-■ Indeed, the important point of our approach is in 
taking into account the uniform part of T, which for quasi-uniform tridiagonal (QUT) matrices is almost the whole 
T, and use the properties of Chebyshev polynomials for reducing the complexity of Eqs. ([5]) and Again, without 
loss of generality we set a — Q and 6 = 1 in what follows. 

The results we present in this paper arc based on the following important statement: the characteristic polynomial 
of QUT matrices can be expressed in terms of the Chebyshev polynomials of the first and second kind, Tn+i{C) 
and Unii), 

X(20 ^ X,-.m = u{S,)Un{0+mTn+l{S.) , (12) 

where u(^) and t{^) are low-degree polynomials: indeed, their degree cannot be larger than i—n and n— 1, respec- 
tively. Their coefficients involve the nonuniform matrix elements and generally they can be easily calculated by means 
of Eqs. Q. 

In order to prove the above general statement we start from the characteristic polynomial of the uniform tridiagonal 
submatrix T^-.v and calculate the characteristic polynomial of larger submatrices by means of Eq. (|4ap : 

X„:„(2a=Z^n(0 

X„^^^, (2^ = (2^ - X„..+, (2e) - 6^+1 Ur^iO 

x„,,(20 =Po(ow«(e) +pmun-m ; (13) 

this holds for some polynomials po{C) — (2^)^^" -I- . . . and Pi(C)i whose coefficients are products of the nonuniform 
matrix elements a^+i, . . . , and b^, . . . ,bi. By further enlarging the matrix T^-i in the upper corner by means of 
Eq. (|4bp we first obtain 

(20 = (2C - a„_i) X.., (20 - bl_, x.^,, (2^) , 

where x,,+i-ei'^0 concerns the QUT matrix Tu+i-.i whose uniform block is (n— l)x (n— 1), so its expression analogous 
to Eq. (fO)) involves Un-i{S,) and Un-2{£.)- Proceeding further one has 

X„_,,(2e) = (2^- a„_2)x„_,,(20 - bl_^x.j20 



X^..M)=P0{mn{O+Pl{mn-m+P2{mn-2{£,) , (14) 
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for some polynomials po(Oi P'^iO^ ^^'^ P'iiC)- This expression allows us to recover Eq. (jl2p . using the identities 

Un-i{0=iU^{0-rn+m . (15a) 

Un-2{£.) = 2ez^„_i(e) -z^„(e) = {2e-i)Un{i) - 2e7;.+i(o , (i5b) 

and identifying 

<0 - Po(e) + CPi(e) + (2^' - l)P2(e) , (16a) 

^(0 = -pi(e)-2^P2(e) ■ (16b) 

The usefulness of expressing x(A) = Xi-.iW the form (|T2|) is evident looking at the analog of Eq. ([T0| for the 
first-kind Chebyshev polynomials, 

7^(cosfc) = cos(nfc) , (17) 

which turns Eq. (jl2p into 

X(2 cos fc) = m(cos k) ^^^^^^nk'^'^ ^ '=o^[("+l)^] ' (1^) 

hence, the secular equation x — ^ can be written 

sin[(n+l)/c - 20fc] = , (19) 

with the angle (pk defined by 

sinfc i(cosfc) 

tan20fc = , ■ 20 

u(cosk) 

Equivalently, the same form of the secular equation can be derived directly by simply rewriting Eq. (jlOp as 
sinfcZ^„(cosfc) = 3{e*'^"+^-''^} and replacing it in Eq. (fTl)) . which turns indeed into 

cj{e»[("+i)fc-20.]| ^0 , (21) 
where 20^ coincides with the phase of the complex number 

zi;, =Po(0 +e-''>i(0 +e"'''p2(0 = kfele""'^'' • (22) 
It is convenient, in order to easily recover the limit of a fully uniform matrix T , to use slightly modified versions 
of Eqs. (HH) and (EOl), namely 

sin[(^+l)A: - 2(^fc] = , (23) 

with shifts (fik defined by 

^„ ^ , 1 t(cosk) sink 

2ipk = (i-n) k - tan-i , \, ■ 24 

u(cos k) 

Hence, the eigenvalues of the QUT matrix, parametrized as A = 2cosfc, with k S [0, vr], can be obtained from the 
equations 

TT j + 2(p kj 

which determine the allowed values of fc. Comparing with Eq. (|lip it appears that the shifts ipk represent the deviation 
from the uniform case, where they vanish. Eq. (j25p can be solved numerically for any j (except for a few j's if there 
are out-of-band eigenvalues, see below). Usually, an iterative computation is fast converging; in the limit of £ S> 1 
even the truncation of (|25p after the first iteration can be very accurate, as it was verified in the cases considered in 
Refs. 

Noteworthy, in the limit of large I Eq. (|25p allows us to obtain a useful analytic expression of the density of states 
Pk defined in the interval fc g [0, tt], namely 

= = ' ^''^ 

by means of which summations over eigenmodes can be transformed into integrals over fc. 



fc = fc, - , , ^ (j = i,...,^), (25) 



Jo 



dkpk {■■■); (27) 



one can also observe that pj, ^ represents the spacing between subsequent allowed values of fc: the deformation from 
the equally-spaced fc's of the uniform case, 7r/(£-f 1), is represented by the correction term with (/j'j,. 
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C. Out-of-band eigenvalues 

The fact of setting A = 2cosfc does not imply that all eigenvalues are included in the band [—2,2]. For a QUT 
matrix this is generally true for the largest part of the spectrum, though a few eigenvalues can emerge over or below 
the band when (the absolute values of) the nonuniform matrix elements are large enough; correspondingly, Eq. (j25p 
cannot be solved for a few values of j, i.e., Eq. ([25]) has less than £ solutions in the interval k G [0, tt]. On the other 
hand, the out-of-band eigenvalues are still described by A = 2cos/c, but with complex values of fc = q + ip; for the 
eigenvalues to be real q must be either or tt, i.e., 

A=±2coshp, (28) 

and p > 0. Correspondingly, one can take the expression for the Chebyshev polynomials when the absolute value of 
the argument is larger than one, 

[/„(±coshp) = (±)»5!!^M^. (29) 

In the large-£ limit, the out-of-band eigenvalues have to be considered separately by adding to the integral ([57]) the 
sum over the out-of-band states. An example of how to deal with such eigenvalues is given in Section llV Al 



D. Eigenvectors 

The boundary elements of the eigenvectors given in Eq. ([5]) can be calculated using the same formalism. Indeed, 
following the construction of the previous subsection we can find the polynomials u (^), t (^), Uj(^), tj{^) such that 

x,,_,(20 = w.(e)z^n(o + i.(e)r„+i(o , (30) 

where the symbols ^ and j clearly refer to the submatrices T2:i and Ti:f_i, respectively. Accordingly, expressing x'W 
as a function of Un and 7^i+i thanks to the relations 

= {n+1) U,,{0 , (31a) 
{l-e) KiO = ^ UniO - [n+l) Tn+m , (31b) 



Eqs. ^ take the form 



, u{£.k)U^{ik)+f{(,k)Tn+l{£.k) 



l2 r,uMk)l^ni^k) +tj{^k)Tn+l{ik) 



<i^k)U„{^k)+t*{ik)Tn+l{^k) ' ^ ' 



where ^k = Afe/2 = cosfc and 



1 

1-e 

n+l 

As the eigenvalues are the solutions of the secular equation 



KiO = u'iO + "(0 + (n+l) m, (33a) 

m)-t'iO-^u{0. (33b) 



- uiCk)Un{^k) + tiCk) r„+i(a) , (34) 
the high-degree polynomials tin{S,k) and Tn+i{S,k) can be removed from (j32p and accordingly 



,2 (6)<(a)-i (a)^6) 



^^'"^<(a)t(a)-i^ (6)^(6) ' ^^^^^ 
2 ^ , ^.fe)tfo)-t.(a).(y (35b) 

K{^k)t{^k)-t*Aik)u{ik) 
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This shows a remarkable result, namely that, although the eigenvector components generally depend on complicated 
high-degree polynomials, for QUT matrices one can express the boundary coefBcients of the eigenvectors in terms of 
ratios of low-degree polynomials. 

Further simplifications can be obtained by replacing again £_k = cos/c. In fact, from Eq. 

„ , ,„ , tu cos k + (t'u — u't) sin^ k 

2(^'fc = (i-n) ^ . ' , (36) 

-f- sm k 

where the argument of u and t is understood, so that the eigenvector elements ((35|) read 



2sin^fc ui^k)ti^k)-fi^k)ui^k ) 



2 _ 2 sin^ k u^{£,k)t{^k) - tj{^k)u{S,k) , . 

^2(^^)+i2(^^) 3in2fc ■ ^^^^^ 

These expressions generalize what was found in Refs. 

As for the remaining elements, note that the recurrence relation ([5]) in the bulk, i.e., for u < v < v, reads 

Ofc,.+i + Oi,.-i = ie'''- + e-'*'-) Ok. , (38) 

whose generic solution is 

Ok. = Qie^'^ak} , (39) 

for any complex number ak independent of which has to be determined by requiring that the 'boundary' relations ([5]), 
i.e., for v = 2, . . . ,u and i> — v, . . . , be satisfied. 



IV. EXAMPLES 



A. Per symmetric two-edge matrix 



As a first example we consider a mirror-symmetric QUT matrix with two non-uniform edges: the uniform block is 
of size n = £—2, so the matrix reads, setting b = 1 and a = 0, 



(X y 
y 1 
1 1 

1 ■•• 



V 



1 

y 

y X, 



(40) 



and, with the notations of the previous section, u = 2 and v = £—1. 

Keeping the notation A = 2^, thanks to the recursion relations (|4]) it holds that 

X,.,M) = {2^~x)U^{0 - y^Un-m 
X,.A2i) - (2e-x)W„_i(0 - y'W„_2(C) 

x,,(2C) = (2e-x)x..(2C)-yV,. 

= {2^~xfUn{S.) - 2{2^-x)y''Un-i{0+y^U^-2{0 . 



(41) 



Accordingly, the secular equation for the in-band eigenvalues is given by (|23p . where the shifts are more easily found 
from Eq. (j22p : indeed, thanks to mirror symmetry, Wk turns out to be a square. 



Wk = {2S,~ x~ y^e"*'')^ = [(2-y^) cos fc - a; -t- iy^ sin k] ^ , 



(42) 
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so that 

^^-'-'^'^ i2-y%cosk-. - ^''^ 

The expression (|^T|) can be rewritten in the form by means of the properties ([T5|) . so with the notation of the 
previous section we identify the coefficients of Eqs. ([T^ and as 

zi(a= [(2VK-2:]'-y*(l-^^) , (44a) 

m = [(2-y2)^ _ ^ (44b) 

"'(e) = w.(e) = (2-y')e-2: , (44c) 

t'iO = tM^y'- (44d) 

Of course, Eq. can be obtained using straightforward trigonometric identities also from and the above 
polynomials. As for the first components of the eigenvectors, they follow from Eq. (|37| : 



^+l-2^'j^ [(2-2/2)cos/fc-x]2 + j/4sin2^ 
Moreover, imposing to the generic solution pop the conditions ([5]) at the corners, one finds 



l„a.e-'fe + (l-y2)e-2'fc 

afc = Oki , (46) 

ysmfc 



so that all the components of the eigenvectors have a fully analytical expression: 



^ sin i^fc — a; sin(j/--l)fc + (1— y^) sin(^— 2)fc ^ ^^^^ 
y sin fc 



for V = 2, 1. Eqs. (|47|) and (j45p . together with (j43p and (l23l) . give a complete solution to the analytical 

diagonalization problem of the matrix . Note that for x — this expression is in agreement with Ref. [IQi] and 
that the term with ip'j^ becomes irrelevant for large n. 

We remark that as long as there are no out-of-band eigenvalues Eq. (1451) is exactly normalized, i.e., J^k ^fci = 1> a- 
sum that in the large-£ limit turns into the integral 

^(^'^)=/ V wo 2, \ 12 a, 4 ■ 2, =1 ■ (48) 

Jo ^ [(2— 2/ ) cos k — xY + sm k 

Eigenvalues A ^ [~2,2] can exist for large x or y. Let us consider the simpler case y — \, with a; > 0, for which 
numerical results are shown in Fig. [1] it appears that when x is raised above 1 two eigenvalues can leave the band. 
Indeed, from Eq. (HTj) one finds the secular equation 

Ut-2xUi-i+x'^Ui-2 = 
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FIG. 2: Eigenvalues of the matrix (|40|) for a; = as a function of y, for matrix sizes ^ = 10 and 50. When y > V2 there can be 
two pairs of opposite out-of-band eigenvalues. 



that by means of the representation ((29l) is easily shown to imply 



X — coshp + sinhp 



tanh 



{£-l)p 



±1 



> 1 



and two out-of-band eigenvalues can indeed exist if x > 1: in the large-^ limit they converge to the same value 
A = X -I- . On the other hand, for a; < 1 all eigenvalues belong to the band. This fact is reflected in the 

z^'^ < 1 for X > 1: indeed, the full normalization requires 



integral because I{x, 1) = 1 for a; < 1, while T{x, 1) 
the contribution from the out-of-band components. 

A similar reasoning applies when y is left to vary while a; = 0, reported in Figl2j though in this case the out-of-band 
eigenvalues occur as two pairs of opposite sign. We find indeed 1(0, y) = 1 for y < V2, while 1(0, y) = (y^-l)-! < l 
for y > v^. It is to be noted that calculating the integral (|48|) in the general case is not trivial; nevertheless, we can 
say that it evaluates to 1 as long as all eigenvalues belong to the band [—2, 2]. 



B. More mirror symmetric elements 



As a second example we consider a mirror-symmetric matrix with more nonuniform elements on the edges, 

/O X \ 

X y 
V 1 
1 1 



1 1 
1 y 
y X 
X 0/ 



Using straightforward algebra we find 



Wk = [2 - - a;^ + (2 - y^) cos 2k + iy"^ sin 2k] ^ , 

uio = [2ei2-y')~x']'-{i-e)4ey\ 
m-^^y' [2ei2-y')-x'] , 

uiO = uM=^[- + ^'(-2 + y^) + 4(2 - 2,/ + y^)e] 

i{0 = U0 = y'[^{2-y')e-^'] , 



(49) 
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and accordingly 



where = 2 



y ■ 



^Ik — '^M — 



fk — 2k — tan ^ 
2 



sin2fc 



_z2 + (2 - 2/2) cos2fc 
x^y'^ sin^ k 



^+l-2V5fc [z2 + (2 - y2) cos 2k] + y4 sin^ 2k 



(50) 
(51) 



C. Non-mirror-symmetric matrix 



In order to connect our formalism with the results of Ref. [l| let us consider the following non-mirror-symmetric 
matrix 



/x y 
y 1 
1 1 
1 1 



1 

1 1 
1 zJ 



(52) 



where £ = n + 2. We find 

m^2^-x-yh 



ui0^ix-20iz~0 + y'iz^- 



and in particular 



tan 2(j)k 



{x + y^z — 2 cos k) sin k 



(x — 2 cos fc)(z — cos fc) -t- y2(2: cos fc — 1) ' 



(53) 



from which the spectral decomposition follows. In fact, it can be shown that, once Oki is calculated with Eg. (1371) . 
the remaining eigenvectors are given by (H7|) . except for the Ith. one that follows from Eq. ([5]). Ec^. ([5^ extends the 
results of Ref. [l|: for example when a; = 0, y = 1, and z = — 1 we find 2</)fe 



-|fc and 



' ~ 21+1' 

recovering Theorem 1 of Ref. 1]. With the proper parametrization it can be shown that the other theorems of Ref. [l| 
concerning symmetric tridiagonal matrices follow as well. 



V. CONCLUSIONS 



We have introduced a technique for the analytical diagonalization of large quasi-uniform tridiagonal (QUT) matrices. 
The quasi-uniformity has been exploited to show that almost all eigenvalues belong to the same band of those of the 
fully uniform matrix, A = 2cosfc, with k e [0,7r], and that their distribution is a deformation of the equally spaced 
/c's of the uniform case, characterized by shifts (pk, as Eqs. and (|25p show. The first components Oki of the 
normalized eigenvectors are written in terms of ratios of low-degree polynomials p7p that can be easily calculated 
from the non-uniform part of the QUT matrix, while the other components are constructed recursively from Oki using 
Eq. ([5]); exploiting the uniform-bulk property, i.e., using Eq. ([5^ . all components can be expressed as Oki times a 
combination of Chebyshev polynomials, as shown in a particular example by Eq. ([7]). 

In the case of a large QUT matrix the eigenvalues can be described in terms of a modified density of states within 
the band of the corresponding uniform matrix. A limited number of out-of-band eigenvalues can exist and have to be 
accounted for separately as discussed in section IIII CI and exemplified in section IIV Al 
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